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The objective of this paper is to define empirical parameters (or closure models) for an 
initial thermal non-equilibrium porous-media model for use in Computational Fluid 
Dynamics (CFD) codes for simulation of Stirling regenerators. The two CFD codes 
currently being used at Glenn Research Center (GRC) for Stirling engine modeling are 
Fluent and CFD-ACE. The porous-media models available in each of these codes are 
equilibrium models, which assume that the solid matrix and the lluid are in thermal 
equilibrium at each spatial location within the porous medium. This is believed to be a poor 
assumption for the oscillating-flow environment within Stirling regenerators; Stirling l-D 
regenerator models, used in Stirling design, use non-equilibrium regenerator models and 
suggest regenerator matrix and gas average temperatures can differ by several degrees at a 
given axial location and time during the cycle. A NASA regenerator research grant has been 
providing experimental and computational results to support definition of various empirical 
coefficients needed in defining a non-equilibrium, macroscopic, porous-media model (i.e., to 
define ‘■‘closure” relations). The grant effort is being led by Cleveland State University, with 
subcontractor assistance from the University of Minnesota, Gedeon Associates, and 
Sunpower, Inc. Friction-factor and heat-transfer correlations based on data taken with the 
NASA/Sunpower oscillating-flow test rig also provide experimentally based correlations that 
are useful in defining parameters for the porous-media model; these correlations are 
documented in Gedeon Associates’ Sage Stirling-Code Manuals. These sources of 
experimentally based information were used to define the following terms and parameters 
needed in the non-equilibrium porous-media model: hydrodynamic dispersion, permeability, 
inertial coefficient, fluid effective thermal conductivity (including thermal dispersion and 
estimate of tortuosity effects), and fluid-solid heat transfer coefficient. Solid effective 
thermal conductivity (including the effect of tortuosity) was also estimated. Determination 
of the porous-media model parameters was based on planned use in a CFD model of 
Infinia's Stirling Technology Demonstration Convertor (TDC), which uses a random-fiber 
regenerator matrix. The non-equilibrium porous-media model presented is considered to be 
an initial, or “draft,” model for possible incorporation in commercial CFD codes, with the 
expectation that the empirical parameters will likely need to be updated once resulting 
Stirling CFD model regenerator and engine results have been analyzed. The emphasis of 
the paper is on use of available data to define empirical parameters (and closure models) 
needed in a thermal non-equilibrium porous-media model for Stirling regenerator 
simulation. Such a model has not yet been implemented by the authors or their associates. 
However, it is anticipated that a thermal non-equilibrium model such as that presented here, 
when incorporated in the CFD codes, will improve our ability to accurately model Stirling 
regenerators with CFD relative to current thermal-equilibrium porous-media models. 
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Schematic of Infinia’s TDC 55 We Stirling Engine 
Pistons have Clearance Seals=>No Contact with Cylinders 

Uses Helium “Working” Gas 

Heater 






GRC Temperature Contours: Generated with Dyson’s Fluent 
2-D Axisymmetric CFD Model of Infinia’s TDC Stirling Engine (Using 
Thermal-Equilibrium Porous-Media Model for Regenerator Simulation) 
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Currently Used Regenerator Random Fiber Material 
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Evidence that Gas/Solid Thermal Equilibrium Assumption 

is not Valid 
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Regenerator Research Test Section at UMN 
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UMN Wire Screen Regenerator 

(Considered to Yield an Approximation of Random Fiber) 

Regenerator^ 
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UMN Test Data: Temperature Difference between Gas/Solid in Regenerator 
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CSU CFD “Micro-” and “Macro-” Calculations & UMN Test Data 
(Micro-Calculations made using REV of actual geometry) 
(Macro-Calculations made with CFD-ACE Thermal Equilibrium Porous Media Model) 
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Regenerator Gas and Solid Temperatures at Hot-End, Mid-Point and Cold- 
End as Predicted by 1-D Sage Code — for Infinia’s TDC 55 We Engine 
(Using 1-D Thermal Non-Equilibrium Regenerator Model) 
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Non-Equilibrium Porous-Media Conservation Equations 

Continuity: j 8 (pY ^ 1 „.L\//..\l_n i 
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Determination of Permeability and Inertial Coefficient 

(1) U. of Minn. Large-Scale Wire Screen Measurements 
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Fluid & Solid Effective Thermal Conductivity Tensors 
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Estimates of Fluid-Stagnant and Solid-Effective Thermal Conductivities 
--for 90% porosity wire screen or random fiber matrix — #1 
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Estimates of Fluid-Stagnant and Solid-Effective Thermal Conductivities 
--for 90% porosity wire screen or random fiber matrix — #2 

• For the axial direction use a modification of the series model 

• A lumped effective solid + fluid effective thermal conductivity (not 
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Heat Transfer Coefficients Between Fluid & Solid Matrix 
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Hydrodynamic Dispersion Term ( Possibly Negligible 
Effect for this Application — This Needs to be Checked) 

• Hydrodynamic dispersion term-fluid momentum equation: — (uu)=uu 
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